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Abstract 

This paper is concerned with the kinetic model of Othmer-Dunbar-Alt for bacterial motion. 
Following a previous work, we apply the dispersion and Strichartz estimates to prove global exis- 
tence under several borderline growth assumptions on the turning kernel. In particular we study 
the kinetic model with internal variables taking into account the complex molecular network inside 
the cell. 
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1 Introduction and results 

In biology, several key processes of cellular spatial organisation arc driven by chemotaxis. The effective 
mechanism by which individual cells undergo directed motion varies among organisms. We are par- 
ticularly interested here in bacterial migration, characterized by the smallncss of the cells, and their 
ability to swim up to several orders of magnitude in the attractant concentration. Several models, 
depending on the level of description, have been developed mathematically for the collective motion of 
cells [■-)?>, '4')]. Among them the kinetic model due to Othmer, Dunbar and Alt (ODA) [1. 'M], describes 
a population of bacteria in motion {e.g. E. Coli or B. Subtilis) [13] in a field of chemoattractant (a 
process called chemokinesis) . These small cells are not capable of measuring any head-to-tail gradient 
of the chemical concentration, and to choose directly some preferred direction of motion towards high 
concentrated regions. Therefore they develop an indirect strategy to select favourable areas, by detect- 
ing a sort of time derivative in the concentration along their pathways, and to react accordingly [27]. 
In fact they undergo a jump process where free movements (runs) are punctuated by reorientation 
phenomena (tumbles) [>;n]. For instance it is known that E. Coli increases the time spent in running 
along a favourable direction [27, 4, 13]. 

This jump process can be described by two different informations. First cells switch the rotation 
of their flagella, from counter-clockwise CCW (free runs) to clockwise CW (reorientation, or tumbling 
phase), and conversely. This decision is the result of a complex chain of reactions inside the cells, driven 
by the external concentration of the chemoattractant [15, 37, 38]. Then cells select a new direction. 
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VC is grateful to University of Edinburgh for the kind hospitality during a one week visit. 
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Although we expect large organisms (like the slime mold amoebae D. discoideum) to choose directly 
a favourable direction, bacteria are unable to do so, and they randomly choose a new direction of 
motion. Actually some directional persistence may influence this selection, privileging some angles 
better than others. However we will not consider inertia here for simplicity. 

From the molecular point of view, the frequency of tumbling events is driven by a regulatory 
protein network made of the membrane receptor complex (MCP), the switch complex located at the 
flagella motor, and six main proteins in between (namely CheA, CheW, CheY, ChcZ. ChcB and 
CheR - more arc involved in B. Subtilis, but the whole picture is similar [15]). This regulatory 
network exhibits a remarkable excitation/adaptation process [4, ■'-!(>, .37] . When attractant concentration 
increases suddenly the tumbling frequency decreases in a short time scale (excitation), but increases 
back to the basal activity after a while (adaptation) . This allows bacteria to follow favorable pathways 
over several orders of the concentration magnitude. Note that a similar adaptation process is involved 
in bigger organisms like D. discoideum [20, 32]. Realistic models have been proposed based on the 
complete regulatory network [17, 37], as well as toy models capturing the key behavior (basically made 
of a two species relaxing ODE system [14]). Note that this network is also known to select positive 
perturbations of the chemoattractant concentration only [4], and to be highly sensitive to very low 
changes in the chemoattractant concentration [ *'■]. 

As a drift-diffusion limit of the ODA kinetic model, one recovers the so-called Kellcr-Segcl model [18, 
7, 6] , where diffusion and chemoscnsitivity coefficients can be derived from the mesoscopic description. 
The Keller-Segel model exhibits a remarkable dichotomy where cells aggregate if they are sufficient 
enough, and disperse if not [21]. Particularly in the two dimensional case, the total mass of cells is the 
key parameter which selects between these phenomena (respectively global existence versus blow-up in 
finite time). This simple alternative is depicted in the whole space in [2]. In the three dimensional 
case however, the relevant quantity ensuring global existence is rather the L^^^ norm of the initial cell 
density [9]. Therefore it is of interest to ask the question of global existence at the mesoscopic level. 
As far as we know, no blow-up phenomenon has been found in the ODA kinetic model. 

The goals of this paper are the following two. First we investigate global existence theory for several 
kinetic models depending on the growth of the reorientation kernel with respect to the chemical. In a 
previous work we succesfully applied dispersion and Strichartz estimates to kinetic models including 
delocalization effects [:]], that can be either a time delay effect due to intracellular dynamics, or 
some measurement at the tip of a cell protrusion. Those techniques are applied here to a class of 
assumptions where the reorientation kernel is actually independent of the (inner and outer) velocities. 
Those assumptions are very rough from the biological point of view, but they aim to determine the 
critical growth of the turning kernel ensuring global existence. On the other hand we apply those 
ideas to a more realistic kinetic model including internal molecular variables, improving the results 
of [12]. We present general assumptions for global existence that can be satisfied by the two species 
excitative/adptative ODE system, or more generally by a complex network. 

We consider the following ODA kinetic model for bacterial chemotaxis: 

dtf + v-V,f= f T[S]it,x,v,v')fit,x,v')dv' 
J v'ev 

-/ T[S]{t,x,v',v)f{t,x,v)dv' , t>Q,x£W^ (la) 

-~AS + S = pit, x) = [ f{t,x,v)dv , (lb) 

associated with the initial condition f{0,x,v) = fo{x,v). The space density of cells is denoted by 
p{t, x). We assume in this paper that the space dimension is d = 2 or d = 3. We assume as usual that 
the set y e of admissible cell velocities is bounded. The free transport operator dtf + v ■ V xf 
describes the free runs of the bacteria which have velocity v. On the other hand, the scattering 
operator in the right hand side of (la) expresses the reorientation process (tumbling) occuring during 
the bacterial migration towards regions of high concentration in chemoattractant S. 

Partial review of plausible reorientation mechanisms. 

We review below the assumptions existing in the literature concerning the reorientation kernel, in order 
to motivate the forthcoming work. 
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Delocalization effects. 



In a previous article [3], we considered mild assumptions of the type: 

<T[S]it,x,v,v') < C (l + S{t, X - ev') + \V S\{t, X - ev')j , (2) 

or, 

< T[S]{t,x,v,v') < c(l + S{t,x + ev) + \VS\{t,x + ev) + \D^S\{t,x + ev)y (3) 

Those assumptions were studied for example in [7, 22] in two or three dimensions of space. 

Under assumption (2), the bacteria take the decision to reorient with the probability X[S] = 
J T[S]{t, x,v' ,v) dv' , and then choose a new direction randomly. Therefore the turning frequency 
increases once cells have entered a favourable area, say (where some delay effect due to internal dy- 
namics is expressed by the space shifting —ev'; the concentration measurement is performed at position 
x — ev' by the cell with velocity w', turning at position x). Intuitively, the cells increase the turning 
frequency to be confined in highly concentrated areas. 

The hypothesis (3) is even more intuitive: cells, when they decide to turn (due to a complex 
averaging over the surrounding area within a radius ~ e), simply choose a better new direction v with 
higher probability. This anticipation measurement can be the result of sending protrusions in the 
surrounding, or considering that the cells have some finite radius with receptors located all over the 
membrane (see also [19] for a similar interpretation at the parabolic level - volume effects have also 
been considered at the kinetic level in [8]). However this interpretation is hardly relevant for bacteria 
which are small cells, unable to feel gradients and to send protrusions. 

Remark 1. The gradient in assumption (2) has to be motivated because we highlight that bacteria 
cannot feel gradients of chemical concentration. As a matter of fact, from a homogeneity viewpoint, 
VS" has the same weight as the time derivative dtS. Therefore we can replace indeed (2) by the 
assumption 

< T[S]{t,x,v,v') < c(l + S{t,x - ev') + \dtS\{t,x - ev')'^ , 

which makes sense biologically (although a more realistic assumption is expressed for example in [10], 
see below (5)j. To see that VS* and dtS do have the same homogeneity, observe that 

dtS = G*dtp = -G*W ■ j , 

where G is the Bessel potential, and the flux j is given by 

j{t,x) = / vf(t,x,v) dv , \jit,x)\ < (max|u|)p(t,a;) . 



As a consequence, in the three dimensional case we have 
\dtS\ 



The dispersion lemma turned out to be a powerful tool for dealing with those assumptions (even the 
second derivative of S can be added in (3)). It turns out that putting together those two hypotheses 
(2) and (3) is a much harder task (due to the fact that we loose the benefit of the decay term in the 
balance law along the estimates). Some progress in this direction was recently made in [7, 22, 3] but 
the whole picture is not clear so far. For example in [3] it was shown that in d = 3 dimensions we have 
global existence of weak solutions if 

< T[S]{t,x,v,v') < l + S{t,x + sv) + Sit,x~ev') + \VS{t,x + sv)\ + \V S{t, x ~ ev')\ , (4) 

provided that the initial data are small in the critical space L'^/^. If (4) is strengthened by dropping 
the last term, then a global existence result was established without a smallness assumption on the 
initial data. The proofs use the dispersion and Strichartz estimates of [•">] and rely on the delocalization 
effects induced hy x + ev and x — ev' . 

Interestingly the fact that some directed motion emerges from turning kernels which resemble to 
assumptions (2), (3) or (4) ~ as pointed out by the diffusion hmit - seems to involve a completely 
different mechanism from the following commonly described behaviour in E. coli. 
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Persistence of motion in the good directions. 

As opposed to the previous set of hypothesis, it is commonly accepted that bacteria increase the time 
spent in running in a favourable direction [38, VS]. That is, the turning kernel is expected to decrease 
as the chemical concentration increases along the cell's trajectory, like 

T[S]{v,v') = To + ij{St+v' -WS) , (5) 

where t/; is nonnegative and decreasing, and St + v' ■ VS* denotes the directional derivative along the 
free run before turning (see [12], [10] where this hypothesis is injected in a model for D. discoideum 
self-organization, and its drift-diffusion limit is derived). One may think of tp to be: ipiv) = if 77 > 
and ip{ri) = 1 if 77 < for instance. Actually, in [12] the authors explicit two behaviour caricatures, 
where cells might "perfectly avoid going in wrong directions", or "perfectly follow good directions". 
The latter is stressed out there and leads the system to regular solutions, intuitively, whereas the 
former might develop singularities where cells aggregate. 

The above mechanism is also part of more complex models including internal variables (which is 
reviewed and analysed further below) . In fact some molecular concentration denoted by y (standing for 
the phosphorylated CheY-P) which induces a tumbling behavior, is actually reduced under attractant 
binding to the membrane receptor (excitation phase). The chemical chain of reactions is in fact 
inhibited under activation of the membrane complex receptor. On the contrary, expression of a repellent 
activates this internal network, favouring tumbling. Global existence theory for such a class of models 
has been discussed in [12] for the one dimensional case. 



Internal dynamics 

Complex models of bacterial motility include a cascade of chemical reactions. This chain of activa- 
tor/inhibitor reactions links the evaluation of the chemical concentration by the membrane receptors 
to the rotational switch of the flagella, inducing or inhibiting the tumbling phase. Several works 
propose a chemical network describing this complexity [17, 37]. In particular, the global short term 
excitation/mid term adaptation is crucial for the cells to crawl up across levels of magnitude of the 
chemical concentration. Caricatures of such an excitation/adaptation process are depicted in [13, 10] 
for instance. However we will keep in this paper the necessary abstract level required for our purpose 
( for an illustrative example, see section 5). 

In the following, y E M™ denotes the whole internal state of the cells, which can correspond to 
huge data of molecular concentrations in the chemical network (in fact m = 2 in the caricatural 
excitating/adaptating system). In accordance with previous notations, p{t,x,v,y) denotes the cell 
density at position x, velocity v, and with internal state y. As before, f{t, x,v) = J p{t, x, v, y) dy is the 
cell density in position x velocity space. On the other hand we introduce ^i{t, x, y) = J^p{t, x, v, y) dv, 
and as usual p{t,x) = yp{t,x,v,y) dvdy. The chemical potential is given by a mean-field equation 
—AS* + S = p{t, x). But this could be extended to a more realistic influence of the internal state on 
the chemical secretion (as it is in [lo]) 

-AS + S= / Lu{y)n(t,x,y) dy , 

under suitable assumptions on the weight uj. The dynamic inside an individual cell is driven by an 
ODE system representing the protein network in an abstract way: 

^ = G(y,5(t,x)) , yeR™. 

The cell master equation describing the run and tumble processes, and the chemical potential equation 
are respectively: 

dtp + v- V^p + Vy ■ (G{y, S)p) = / T{t, x, v, v', y)p{t, x, v', y)dv' 

V / Jv'eV 

- T{t,x,v',v,y)p{t,x,v,y)dv' , (6a) 

Jv'ev 

-AS + S = p , (6b) 
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The turning kernel T can be decomposed in this context as product between a turning frequency 
X[y], depending on the internal state only, and a reorientation K[v,v') which may describe some 
persistence in the choice of a new direction with respect to the old one. Without loss of generality 
here we assume that K(v,v') is constant and rcnormalized as being K[v,v') = 

It is worth noticing that this realistic kinetic model may contain enormous informations on the 
microscopic cell biology, and links different scales of description, because we eventually end up with a 
cell population p{t,x). 

As a partial conclusion, we observe that several scenarios with different underlying kinds of hy- 
potheses, drive the system to positive chemotaxis (at least considering the formal drift-diffusion limit 
of those). 



Statement of the main results 

In this paper we investigate the critical growth of the turning kernel in terms of space norms of the 
chemical which ensures the global existence for the kinetic model. In particular we consider control 
on the turning kernel without any dependence upon the velocity variables, that is with some abuse of 
notations: 

< T[S]{t,x,v,v') < T[S]{t,x) , 

under suitable conditions on the growth of T[S]. We exhibit examples in 2D and 3D, restricting 
ourselves to some norms of the chemical (and not of its gradient for instance) for which our method 
appears to be borderline. In particular, it is natural to ask (see Section 3 of [7] and the concluding 
remarks in [•!]) whether global existence can be established under a hypothesis of the form 

<T[S]it,x,v,v') <c{l + \\S{t,-)\\l^^^,^) , (7) 

where a > 0. 



Exponential growth in dimension 2 

Consider first the case of dimension d = 2. It is easy to see using the methods of [7] that wc have global 
existence for any exponent a > within (7). In analogy with global existence results for nonlinear wave 
or Schrodingcr equations [23, 24, 29, 30] we can ask whether the turning kernel can grow exponentially: 

< T[S]it, X, V, v') <c(l + cxp [||^(t, Oll^oodi^)] ) . (8) 

We will show that this is actually possible: if < /? < 1 then we have global existence for large data; 
ii (3 = \ we have global existence for initial data of small mass. Our proof requires M < tt, but we 
don't know if this bound is optimal. Also, wc don't know if we may have blow-up for large M or for 
exponents (3 > 1. 

We shall prove the following 

Theorem 1.1. Consider the system (1) in d = 2 dimensions under hypothesis (8) and let 1 < p < 2. 

Assume < jS < 1. If (3 = \ assume also that M < tt, where M = Jy fo(x,v)dv is the total mass of 
cells. Then if fo € L\L^ C\ L\ ^ then (1) has a global weak solution f with f{t) G L^L\, H L\ ^. 



Almost L°° growth in 3D 

Naturally, from the global existence point of view we address the question of a Hi^lloo growth of the 
turning kernel in the case of d = 3 dimensions: T[S\ < C{l -\- |j5'||oo)- Actually it cannot be handled 
using our method in three dimensions so far. Even in the simpler case T[S] < C(l -I- S(t,x)) our 
dispersion method fails. Puzzling enough, if T[S] = C(l + S{t, x)) then a very simple symmetrization 
trick does perfectly the job (see section 2). 

It was noticed in [7] that if a < 1 in (7) then we have global existence (a sketch of the proof will 
be given in Section 4.3). The case a = 1 remains open. In this direction we will use the methods of 
[3] to show that we have global existence under the assumption 

0<T[S]it,x,v,v')<c[l + \\Sit,-)\\l.^^,^) , (9) 
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where < a < and r can be arbitrarily large. Notice that ^^-^ ^ 1+ as r ^ oo, which is coherent 
with the above obstruction. More precisely we shall prove the following 

Theorem 1.2. Let d — 3 and 1 < p < 3/2. Suppose that that the turning kernel T[S] satisfies 
hypothesis (9) for some r and a verifying: if 1 < r < 3, a can be any positive number, whereas in 
case of 3 < r < oo, < a < Then if fo G L^L^ H „ then the kinetic model (1) has a global 

weak solution with f(t) G L^Ll D Ll. ^. 

If we assume that the turning kernel satisfies (7) with a = 1, we can use the Strichartz estimates 
of [5] to show global existence, provided that the critical norm ||/o 11^3/2^^6 is small. 

Theorem 1.3. Let d ~ 3 and assume that the turning kernel satisfies 

< T[S]{t,x,v,v') <c[l + \\Sit, •)|li=o(R3)) . 

Assume also that fo £ L^ (K^ ^) H i^/^ (R^ j,) and that |l/o|1^3/2^j{6 is sufficiently small. Then (1) 
has a global weak solution. 

Internal dynamics. 

We shall prove the following theorem for global existence in three dimensions of space. 

Theorem 1.4. Let d = 3. Assume that the turning kernel has the form T ~ A[?/] x /^(w, w') where K 
is uniformly bounded, and A grows at most linearly: X[y] < C(l + On the other hand, assume 

that G has a (sub)critical growth with respect to y and S: there exists < a < 1 such that 

\G\{y,S)<c{l + \y\ + S^) . 

Then there exists an exponent 1 < p < 3/2 such that the system (6) admits globally existing solutions 
with p e L'P.LlLy. 

2 The dispersion lemma applied to kinetic chemotaxis, and 
the symmetrization trick 

In this section we present a direct application of the dispersion lemma [ri] to system (1). As a conse- 
quence we are led to the following question, which is decoupled from (1): 

Investigate the critical norm for the turning kernel ensuring the bound 

o<r[5]<c(i + ||p(t)||Lp) , 

for p < d' in dimension d. 
The rest of this paper will be devoted to this question of critical growth. 

Lemma 2.1. Assume the turning kernel can be controlled without any dependence on the velocity 
variables v nor v' : 

< T[S]{t,x,v,v') < T[S]{t,x). 
Then, applying the dispersion estimate, we get the following for p Cz [l,d'): 

IIpWIIlp < \\fo{x~tv,v)\\L.Li+\V\^/P f (t-s)-^ f T[S]{s,x)p{s,x) dxds , 

Js=a Jx 

where A = d/p' . 

Observe that the condition d < p' is crucially required here to ensure further the time integrability 
of the right-hand-side. 
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Proof. As usual we represent the solution of (1) as 

f{t,x,v)<fo{x~tv,v)+ / T[S]{s, X — (t — s)v)p{s, X — {t — s)v)ds 

Jo 

Using dispersion we get immediately 

\\f{t,x,v)\\L-^L}<\\fo{x-tv,v)\\i^.Ll+ / T. TUT^ T[S]{s,x)p{s,x) 



(t- s)'i(i-i/p) 



ds 



□ 



As an observation, we state also a second lemma, which is interesting in its own right, but which 
will not be used in the sequel. Following [.34], it claims that a kernel which is symmetric with respect 
to V and v' ensures global existence. It is relevant from the mathematical point of view because we 
consider bounds that do not depend on v and v' . It is biologically irrelevant however in the case of a 
purely symmetric kernel because no directed motion emerges in the drift-diffusion limit [7] . 

Lemma 2.2. Consider the scattering equation, 

dtf + v- V,/ = / (Kit, X, V, v')fit, X, v') - Kit, X, v' , v)f{t, x, v)) dv' . (10) 
Jv 

and assume that K is symmetric w.r.t. v and v' , i.e. 

K{t, X, V, v') = K{t, X, v', v) > 0. (11) 
Then all LP LP— -norms of the density f (1 < p < oo) are uniformly estimated like \\f{t)\\j^p < ||/o||^p . 
Proof. First rewrite (10) using the symmetry property. It becomes 

dtf + v- V^f = / Kit, X, V, v') if it, X, v') - fit, X, v)) dv' . (12) 
Jv 

Next multiply (12) hy pfP-'^it,x,v) to get 

dtfP + v-W^fP =p J fP-\t,x,v)Kit,x,v,v')ifit,x,v')- fit,x,v))dv' 
Integrate with respect to x and v to get 

fPdvdx =P jjj F^\t^ X, v)Kit, X, V, v') if it, X, v') ~ fit, X, v)) dv'dvdx . 
We can symmetrize the latter expression to obtain eventually 



d 
di 



^JlFdvdx 



Kiv, v') {fP-\t, X, v) - fP-\t, X, v')) if it, X, v) ~ fit, X, v')) dv'dvdx 



Since / > we have 

{fP~\t,x,v) - fP-\t,x,v')) ifit,x,v) - fit,x,v')) > 
because these two factors always have the same sign. It follows that 



d_ 
di 



fPdvdx < . 

□ 
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3 Exponential growth in L°° in dimension 2 

In this Section we prove Theorem 1.1. Working as in the proof of Trudinger's inequahty we expand 
the exponential into a power series and use Young's inequality as in [7, 'A] to estimate each term. The 
dispersion method is then used as in [3] through (2.1). Throughout these processes we keep track of 
the growth of the various constants in order to make sure that the resulting series converges. A similar 
approach has been used in [23, 24, 29, 30] to study nonlinear wave and Schrodinger equations. We will 
need the following two Lemmas. 

Lemma 3.1. Let G{x) ~ e~^^~ e~T^ ^ . There exists a positive constant A such that 

Gix)<A+^\\og\x\\ , |a;|<l. (13) 
Proof. Fix X with < 1. Write G{x) = Gi{x) + G2{x) + Gsix) where 

47r Jo s 
^ , , If' ^^ds 



For Gi use e ^^^'^ < 1 and then change variables s ^ t, where s = \x\'^t, to get Gi{x) < Jq e ^t^—: 
Ai. For Ga we have G2ix) < ^ J^^^^ ^ = ^^1^- For G3 use e"'^!^!'/'* < 1 to get Ga^x) < 
^J^e-^ds =: A2. 



4tT 7|^|2 s ' 

G3(.x) = -y^ e = e 



□ 



Remark 2. In fact the exact asymptotics of G near the origin is: 

G(a;) = --^log|a;|+7+-^log2 + o(l) , 

ZTT ZTT 

where 7 is the Euler constant. 

Lemma 3.2. For x > define T[x) = t^~^e~^dt. Then (Stirling's formula) 

n! = r(?i+ 1) - y2^(^-j (n^+00), (14) 
T{x + l) ^ V2^(^^y (x^+oo). (15) 



Moreover, for all P > 0, x > 0, 



Proof. (14) and (15) are well known. For (16) we have 

poo poo pOO 

T{f3 + 1) = / t^e-*di > / tf^e-'dt > x'^ / e^'dt = x^e"^. 

Jo Jx Jx 

□ 

Proof of Theorem 1.1. Recall from section 2 that a control of the turning kernellike T [5] < G(||p(t)||Lp) 
is sufficient to guarantee global existence. The rest of this section is devoted to the proof of this esti- 
mate. 

Pick 1 < p < 2 and set ^ ~ p' > 2. In case of /3 = 1, assume in addition that jj < jj. 
Write S = S''°"f + S"'""'* where 

= (l|,|>iG(x)) * p , and 5^''°'^* = (l|,|<iG(x)) * p . 
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Since < /3 < 1 we have 



II^IIloc, < (ll-s'^ii^ + ||s"''°'*IL^y < \\s^°''s\\1^ + \\s 



'short 1 1 ^ 

\\l° 



IL°= II \\L°^I — II IlL" 

where we have used the fact that [x + < + for x,y > and < /3 < 1. Therefore 



exp 



For 5'°"^ we have 



{||5(i,-)ll^=o}<cxp{||5'°"^(t,.)||5!o.}-exp{||5^'""-*(t,.)||?=.} 



where c is a positive constant (depending on /?), therefore 

exp[\\S{t, } < e-'^' exp { \\S^^°^\t, } 



1 1 oshort ( I 



1 + ^- 



For 5"''°'^* we have 



IIpI 



L f J- 



For all 7 > 1 we have , < — ^ = p therefore 

||5^'°'^*(i,-)|Lo. < ||G(x)l|,|<i||^,, Mt,-)\\l~J \\p{t,-)\\l 



therefore 
Consequently 

exp 



||G(a;)l|,|<i||^„ Afi-7||p(t,.)||i,, 



'jshort 1 1 



{\\s{t,-)\\1^} 



< e 



.Ml' 



f;i||G(x)l|,|<i||ll.M^-'' 



LP 



(17) 



We need to guarantee that the series in the above right-hand-side converges. Using (13) we have: 



||G(x)l|,|<i||^„ < 



+ — lloglxl 



< Att H j|log \x 

27r 



Lf'j(|2;|<l) ' 



and also 



||log|x| 



'L^'^{\x\<l) 



1 \ 1/w 

2t: I { — logr )^'' rdr 







< ( 27r / s^^e-^'ds 



< 27r 



^0 r(w+i) 



■e "ds ) by (16) 



As a consequence 



(27r)" (r(MJ + l))" 



|G'(^)lbl<i|L« < + ^(2^)" (r(Aij + l))" . 



(18) 
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Then the infinite sum in (17) can be estimated by 

j/3 



^ - Mtt^ + — (27r)^^ (r(/ij + 1))^^ j AP^. (19) 



We'U show that for /3 < 1 the series converges for any mass M, and that for /? = 1 it converges thanks 
to the restriction 4r^ < 1. Using the root test we have 



1 / 1 1 1 ^V^ rA' 

- Mtt" + — - (27r)" (r(^j + 1))" j AP^ \ 
= 5- M^T" H (27r)" (r(^j + 1))" M'^ 

(j!)7 V 2^ ; 

< Ia^ttTU + ("-L) (27r)^ (r(Mj + 1))" I Ad", 



1 |3 

We have — ^-A'^^ttm ^ 0, therefore it remains to examine the hmit of 



1 fl-" 



(j!)7 V2^y 



(27r)" (r(^j + 1))" M^. 



(20) 



From (14) we have j! ^ V^tFJ (i)"' therefore (j!)^ ^ {2TTj)^i i ^ i- From (15) we have r(/xj + 1) 
y/2^j{ifY^ therefore 



(r(^j + l))" - (27r/xj)^" 



Therefore 

/3 Ciii^'^ 



■\ /5 / ■\P 



(20) ^r^v m:^ . if /5<i 

^ ^ I , if /3 = 1 

The hmit is smaher than 1 in aU cases, therefore the series converges. 
Summing up, we obtain 



\Lp 



T[S]{t, x)<c(l + cxp {|15(i, ])<C + C \\p{t, ■% 
Recall that we have choosen p < 2 such that the Lemma 2.1 applies. We end up with 

\\pit,x)\\^,<t-'\\fo{x,v)\\^.LP+C (l + ||p(s,x)||^.)- ' 



where A = 2/p' < 1 so that we can bootstrap. □ 

4 (Almost) growth in dimension 3 

4.1 Almost L°° growth 

In this Section we consider the kinetic model (1) in d = 3 dimensions under hypothesis (9). 

Proof of Theorem 1.2. Ifl<^<3, a>0 and T[S] satisfies (9) then T[S] can be estimated a priori 
in terms of the mass M. Indeed, 

||^(i,-)llL.(E3) < ||G||^.(K3) \W^-)\\LHm ^ 

because G[x) ~ ^ for smaU and G{x) decays exponentiaUy for large x. Therefore T[S]{t, x, v, v') < 
C + CA/P and global existence follows easily. 



10 



Assume now that 3 < r < oo and < a < ;j7^ . Choose p defined by 



1 a(r~3) 1 

- < 



3r 



3 ' 



and define B sucli that 



11 11 

B' 3 r ap' 

Using fractional integration [20] we get for the signal S = G * p < * p[x) (both short and long 
range parts), 



\\S\\l.<C 



<C\\p\\^s<CM^-'t\\p\\l'^ 



Consequently we get the crucial estimate required in Lemma 2.1: 

T[S]{t,x)<C + C\\S\\l.<C + C\\p\\^, , 
where p is smaller than 3/2. We can complete the proof as in Theorem 1.1. 

4.2 L°° growth: global existence for small data 

Proof of Theorem 1.3. We have 

dtf + v-^.J<cj^{l + \\S{t)\\^^.)f{t,x,v') dv' = C(\ + \\S{t)\\L^ )p{t.x). 

To apply the Strichartz estimate [5] we need four parameters q^p, r, a such that 

1 < r < p < oo 

1 1 1 

0< < - 

r p 6 

1 1 
1 < - + - 

r p 

^^3^-1 
q \r p 

2pr 



p + r 



More conditions will be imposed later. We get: 

\\f\\^LUz<\\foh^_+C il + \\Sit)\\^Mt,x) 
= ll/olL. 



Ci\V\) {l + \\S{t)\\^^.)\\p{t,x)\\^^ 



LI 



□ 



(21) 

(22a) 
(22b) 

(22c) 

(22d) 

(22e) 



(23) 



In the sequel we omit the constant part in the growth of the turning kernel for the sake of clarity. 
Assume 

P>1- (24) 



Then p' < 3 therefore, 



\\S{t)h^. <\\G* p{t)\\^^ < \\G\\^,, ||p(i)ILp < G MmL. 



(25) 



because G{x) ~ -j^ for small \x\, and G{x) decays rapidly for large Moreover, since r < p we have 
by interpolation. 
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Therefore 



\\smL^Mt,x)\\ 



< c 



IpWIIlp \w)\\i. 



= 11 \W)\\l. II 



1+^ 



IIpWIIlp = ll/(^,^,^^)llLSLi ^^^^(I^I)II/(*'^-'^)IIl£l: 

l|5(t)IL^ l|p(t,^)llL. < C||/(i,X,«)ll'tL.'^ 



Now 
therefore 

Suppose that 
Then 

I \\S{t)\\^^\\p{t,x)\\^^^ <C||/(i,x,z;) 
and phigging this into (23) we get 



1 



,1 + 7- 



(26) 
(27) 



ll/(i,2;,t.)IU,iPi. < ||/o|Lj_^+C|l/(t,x,t;)||^,^p^. 

If II/oIIlo is smah enough then we can bootstrap. 

We need to verify that there exist (g,p, r, a) satisfying (22), (24) and (26). There are many possible 
choices. For example, if we want initial data /o G ^, with a = § (critical exponent in dimension 3) 
we must choose p and r so that ^ + ^ = |. The complete set of exponents solving these constraints is: 

, = 1 + V2,p=^±|^, r = 3fV2-l 
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where all conditions are fulfilled. 



□ 



4.3 Sublinear growth 

To close this section we give a quick sketch of the observation in [7] that the hypothesis 



o<r[^](t,x,«,«')<c(i + i|5(t,-)|lioo) , 



implies global existence. Fix p and q such that 

all 3 
3 p g 2 

Then we have the following elliptic estimate (see below), 

||5(t,.)IL^=||G*p(t)||^. <C(Af)||p(Olli'f . 
Therefore (again omitting the constant contribution of the turning kernel) 



f{t,x,v)<fo{x~tv,v) + C / \\p{s)\\1pp{s,x- {t- s)v)ds. 

Jo 



(28) 



12 



Take the L^L^ norm and use the dispersion estimate with \ = 3{l/q — l/p) = a to get 



(29) 



where p'a/3 + p' /q' = 1 by definition. 

To prove the eUiptic estimate (28) write 



Then, if p' < 3, 



X|2;|<ii|l , ,1 ||X|2;|>_R| 



\\Sh^ <C\\ph4^^\\L.' +C\\p\\l^\ 
<c{\\p\\l.R^-' + \\p\\L^Rr') 

Choose R so that \\p\\lpR~ ^ = \\p\\]^iR^'^ , i.e. choose 



R 

This gives 



II \ P'/3 



pWlp, 



5 Extension to internal dynamics 

Recall the kinetic model with internal dynamics: 

dtp + v- + • (G{y, S)p) = I X[y]K{v, v')p{t, x, v' , y)dv' 
^ ' Jv'ev 

- >\y\p{t,x,v,y) , (31a) 
-l\S + S = p, (31b) 

Assuming that K is bounded we reduce to K = 1/\V\ without loss of generality. 

This model takes into account the transport along characteristics of the internal cellular dynamics 

^^Giy,Sit,x)) , yeM". 

For E. coli, the regulatory network described by G is made of six main proteins essentially (named 
Che-proteins), and the main events are methylation and phosphorylation. Indeed, in the absence of 
a chemoattractant (basal activity), the phosphorylated protein Che-Y is supposed to diffuse inside 
the cell and to reach the flagella motor complex, enhancing switch between CCW rotation and CW 
rotation, that is tumbling. This transduction pathway is in fact inhibited when the chemoattractant 
(say aspartate) binds a membrane receptor, triggering methylation of the membrane receptor complex, 
and eventually inhibition of the tumbling process. 

This network exhibits a remarkable excitation/adaptation behavior, which is crucial for cell migra- 
tion. For the sake of simplicity, one deals in general with a system of two coupled ODEs which captures 
the same features. This system should be excitable with slow adaptation - there is a single, stable 
equilibrium state, but a perturbation above a small threshold triggers a large excursion in the phase 
plane (see figure 1 {left)) - and possibly one-sided - in case of positive chemotaxis, the cells do not 
respond specifically to a decrease of the chemoattractant concentration [4]. This characterization of 
dynamical systems is very well known in biological modeling, as it is the basis of the FitzHugh-Nagumo 
models [28] for potential activity in axons. Furthermore, it is often associated to the phenomenon of 
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Figure 1: A two coupled ODE system exhibiting short time excitation/ mid time adaptation with one- 
sided selection. The picture is the same as for the FHN model. The perturbation of the equilibrium 
state is enhanced by a displacement of the basal line y2 = h{SO + dS), translating the system up (left) 
or down (right). Here the cubic function is given by q{u) = u{u — l){u — .2), the saturating ligand 
function is given by h[S) ~ 5/(1 + 5), the basal aspartate concentration is 50 = .4 and the system 
reacts to perturbations dS = .1. 



pulse wave propagation {e.g. calcium waves), [25]. In the context of cell migration, it is also involved 
in the slime mold amoebae D. discoideum aggregation process, where the chemoattractant cAMP is 
relayed by the cells [20, 10]. 

To be more concrete, the following set of equations is generally proposed [13] 

^ = —(^('5') - (yi + y-i)) {excitation) , 

dy2 ^^^^ 
= — {h{S) — y-z) {adaptation) . 

Considered to be decoupled from the transport equation, these two internal quantities relax respectively 
to 

hm yi = , lim y2 = h{S) , 

t — >oo t — >oo 

with a slow time scale associated to adaptation provided that ^ Ta. However, this system cannot 
reproduce true excitability with a large gain factor for small perturbations because it is linear with 
respect to the variable In a slightly different context (pulsatory cAMP waves), Dolak and Schmeiser 
considered an even simpler system [10], namely 

yi = {h{S) —2/2)^ {excitation) , 
~ — {J^{S) ~ y2) {adaptation) . 

This particular choice does select responses to one-sided stimuli, but fails for true excitability. We 
suggest to consider the following phenomenological translated slow-fast, FHN type, system. 



^= —{h{S)-q{yi)-y2) {excitation), 

dt Te 

^ = —{h{S) + y{l) - 2/2) {adaptation) , 

dt Ta 



(34) 



where g is a cubic function depicted in figure 1. 

Proof of Theorem I.4. Our next step is to prove global existence under general and manageable as- 
sumptions settled in Theorem 1.4. But let us begin with an important remark on the methodology. 
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Remark 3. To obtain a priori estimates, one possible strategy would be to use the characteristics to 
handle with (31), as it is performed in [I:.:!] in ID. For this purpose, integrating the hyperbolic (31) 
along the backward-in-time auxiliary problem 



X{s)=v, Y{s)^GiY,Sis,X)) , {X{t),Y{t)) = ix,y) , 



gives the estimate 



ds 



/(s, X{s),v, Y{s)) - (Vy ■ g)p < X[YUs, X, Y) 



The difficulty arises at two levels here. First one has to control the Vy • G contribution, and secondly 
one has to perform later on the change of variables z = Y{y). This induces a Jacobian contribution 
dY 



, and one has to control it too. In the sequel we avoid these two difficulties by working on 
dy 

averaged quantities. 

We use a partial representation formula of the solution from the free transport operator dtp-\-v-'S/xP- 
First integrating the equation with respect to y, we obtain 



dtf + v-V,f + 0<-^^J^ X[yMt, X, y) dy 



so that 



/(t, X, v) < foix -tv,v)-\~ — . 

I ^ I Js=0 



A[y]/i(s, x-{t- s)v, y) dyds 



Using the L'P.Ll dispersion Lemma 2.1 we get as usual 



\\pmLP<\\foix~tv,v)\\LU^+ — 

I ^ I Js=0 



X[y]fi{s,x- {t - s)v,y) dy 



ds 



<t-^\\fo{x,v)\\LiLZ 

+ 1^1'^""' / 77^ 1 1 AMM(s,x,y) dxdyds , 

Js=0 ['■ ^ ^) J Jx,y 

where A = 3/p'. We now use the two growth assumptions on A and G: 

AM < C(l + \y\) , \G\{y, S) < C(l + |y| + 5") , < a < 1 

to control the time growth of the average quantity JJ^ ^ IvlfJ'it ~ x, y) dxdy. 
Remark 4. Note that in dimension d = 2 we can handle any nonnegative a. 
We test the master equation (31) against \y\: 



d 
ds 



\y\^j,{s,x,y) dxdy + + 



\v\'^y ■ {G{y, S)ii{s, X, y)) dydx = 



x,y 



therefore, using |G| < C(l + |y| + 5"), 



d 
ds 



\y\^J,{s) dxdy 



y 



G{y,S)n{s,x,y) dydx 



< 



x,y [Vl 

\G\{y,S)fj,{s,x,y) dydx 



ix,y 

<G + C 



\y\li{s,x,y) dydx -\- C / \S{s,x)\°'p{s,x) dx 



(35) 



(36) 
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Remark 5. If we agree to diminish a it is possible to deal with higher exponents in X[y] < C(l + 
For instance we have by Young's inequality: 



d 
ds 



\y[^ ^-{s) dxdy < ^ I I \yy ^\G\{y, S)fi{s,x,y) dydx 

x,y 'J <J x,y 

\yrKs,x,y) dydx 

J J x,y 

+ Cl II {^^\yV + -\S{s,x)^\^i{s,x,y) dydx . 

The same argument follows provided that 07 < 1 . More general combination of exponents could have 
been considered. We have chosen here a simple framework for the sake of clarity. 



We can use the Duhamel formula to represent the inequaUty (36) as 
/ / \y\f^is, X, y) dxdy < Ce'^' + e<^" / / \y\tJLQ{x, y) dxdy 

J -J x,y J J x,y 



+ C / e^^'-''^ I \S{t,x)\°'p{t,x) dxdT 

Jt=0 



Plugging that into (35) gives 



WpWlp < Co{t)+C / -3^ / e^(^-"' / \SiT,xTp{T,x) dxdrds . 

We choose p < 3/2 so that A = 3/p' < 1. Since a < 1 we have 3 < ^ and we can choose p sufficiently 
close to 3/2 so that 3 < p' < f . Then 

S{t,xrp{t,x)dx < ilpll^, = \\ph. . 

^From the mean field chemical equation (31b) — A5' + S = p we deduce the following elliptic estimate. 
We have ap' < 3, and S = G * p where G{x) ~ |^ for small \x\ (short range) and G{x) decreases 
exponentially fast for large \x\ (long range). Thus we obtain 



therefore 



We obtain 



= <llG|lw \\ph^<CM 

S{t,x)°'p{t,x)dx <CM"\\p\\^, . 



\\pmL. < Coit) + C / ^ / ||p(r)||^, drds 

Jo ^ *7 Jr=0 



<Coit) + C J^\\pir)\\^, J ^^—^dsdr, 

Using the boundedness of J^^^ {t-s)>' e'^'-^^^-' ds with respect to r, we conclude thanks to a Gronwall 
estimate. 

□ 
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